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We report on the results of an analysis of XMM-Newton observations of 
MS0735+7421, the galaxy cluster which hosts the most energetic AGN outburst 
^sO [ currently known. The previous Chandra image shows twin giant X-ray cavities 

(~ 200 kpc diameter) filled with radio emission and surrounded by a weak shock 
front. XMM data are consistent with these findings. The total energy in cavities 
^ \ and shock (~ 6 x 10 61 erg) is enough to quench the cooling flow and, since most 

of the energy is deposited outside the cooling region (~ 100 kpc), to heat the 
gas within 1 Mpc by ~ 1/4 keV per particle. The cluster exhibits an upward 
^ ' departure (factor ~2) from the mean L-T relation. The boost in emissivity pro- 

duced by the ICM compression in the bright shells due to the cavity expansion 
may contribute to explain the high luminosity and high central gas mass fraction 
that we measure. The scaled temperature and metallicity profiles are in general 
agreement with those observed in relaxed clusters. Also, the quantities we mea- 
sure are consistent with the observed M-T relation. We conclude that violent 
outbursts such as the one in MS0735+7421 do not cause dramatic instantaneous 
departures from cluster scaling relations (other than the L-T relation). However, 
if they are relatively common they may play a role in creating the global cluster 
properties. 
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cellaneous — X-rays: galaxies: clusters 

1. Introduction 

The current generation of X-ray satellites, Chandra and XMM-Newton, has shown that 
the physics of the intra-cluster medium (ICM) is complex and needs to be regulated by 
additional, non-gravitational processes beyond simple gravity and gas dynamics considered 
in the standard Cold Dark Matter cosmological scenario (White & Rees 1978). In particular, 
our understanding of cooling flow systems has radically changed. Albeit confirming the 
existence of short cooling times, high densities and low temperatures in the cluster cores, 
the arrival of new high-resolution X-ray data has shown the lack of the emission lines expected 
from gas cooling below 1-2 keV and has reduced by about one order of magnitude the new 
spectroscopically-derived mass deposition rates (e.g. David et al. 2001, Johnstone et al. 
2002, Peterson et al. 2003 and references therein). The most plausible solution to this 
so-called "cooling flow problem" is that some form of additional heating which balances 
the cooling must be acting in the ICM. Among the many proposed heating mechanisms^, 
one of the best candidates for supplying the energy is feedback from the Active Galactic 
Nucleus (AGN) hosted by the central galaxy of the cluster (e.g., Rosner & Tucker 1989; 
Tabor & Binney 1993; Churazov et al. 2001; Briiggen & Kaiser 2001; Kaiser & Binney 2003; 
Ruszkowski & Begelman 2002; Brighenti & Mathews 2003, Omma et al. 2004). 

The possibility of AGN heating has recently become the leading idea to solve the cooling 
flow problem thanks to the discovery of X-ray cavities in the ICM on scales often approxi- 
mately coincident with the lobes of extended radio emission (e.g., Hydra A: McNamara et 
al. 2000, David et al. 2001; Perseus: Bohringer et al. 1993, Churazov et al. 2000, Fabian 
et al. 2000; A2052: Blanton et al. 2001, 2003; A2597: McNamara et al. 2001, Pollack et 
al. 2005, Clarke et al. 2005; RBS797: Schindler et al. 2001, Gitti et al. 2006; A4059: Heinz 
et al. 2002). This indicates that the radio lobes fed by AGNs have displaced the X-ray 
emitting gas, creating cavities in the ICM. These outbursts can be recurrent, as it has been 



1 Proposed heating mechanisms include electron thermal conduction from the outer regions of clusters 
(Tucker & Rosner 1983; Voigt et al. 2002; Fabian et al. 2002; Zakamska & Narayan 2003), continuous 
subcluster merging (Markevitch et al. 2001), contribution of the gravitational potential of the cluster core 
(Fabian 2003), feedback from intra-cluster supernovae (Domainko et al. 2004), etc. 
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observed in some cases. For example, RBS797 shows evidence of a restarted AGN activity 
with precessing jets pointing to a different direction with respect to the one of the outer radio 
lobes filling the cavities (Gitti et al. 2006). The heating by AGN is thought to occur through 
the dissipation of the cavity enthalpy and through shocks driven by the bubbles inflated by 
the AGN. Systematic studies of a sample of X-ray cavities show that their enthalpies^ lie 
between 10 55 — 10 61 erg and scales in proportion to the cooling X-ray luminosity and the 
radio power of the host system. In more than a half of the sample the energy input from the 
central radio source is currently sufficient to balance cooling (Birzan et al. 2004, Rafferty et 
al. 2006). The trend between X-ray luminosity and bubble mechanical luminosity, together 
with the existence of short central cooling time, suggests that the AGN is fueled by a cooling 
flow that is itself regulated by feedback from the AGN. The basic idea of this AGN-cooling 
flow scenario is that a self-regulated equilibrium may be achieved, in which the energy input 
from the AGN balances the radiative losses of the ICM over the lifetime of the cluster. 

The recent discovery of giant cavities and associated large-scale shocks in three system 
(MS0735+7241: McNamara et al. 2005, Hercules A: Nulsen et al. 2005a, Hydra A: Nulsen 
et al. 2005b) has shown that AGN outbursts can not only affect the central regions, but 
also have an impact on cluster-wide scales. In particular, the supercavities discovered in 
MS 073 5 +7421 have a diameter of about 200 kpc each and a weak cocoon shock surrounding 
them has been detected. This large-scale outburst is the most powerful known so far: it 
releases upward of 10 61 erg into the ICM, heating the gas beyond the cooling region (Mc- 
Namara et al. 2005). In this paper we want to investigate the significant consequences that 
this new development has for several fundamental problems in astrophysics. One important 
problem addressed by giant cavities is the so-called cluster " pre- heating" , which is mani- 
fested in the steepening of the observed luminosity vs. temperature relation for clusters 
with respect to theoretical predictions that include gravity alone (e.g., Markevitch 1998). 
Some extra non-gravitational energy is required to explain such a steepening (e.g., Wu et al. 
2000, Voit 2005), and one possibility is that it is supplied by AGN outburst. The additional 
non-gravitational heating supplied by AGN could also induce the suppression of the gas 
cooling in massive galaxies required to explain the exponential turnover at the bright end 
of the luminosity function of galaxies (Benson et al. 2003). This would indicate a common 
solution for the two major heating problems associated with the ICM: those of cooling flow 
and galaxy formation. Some cavities live longer than the time it takes to cross their own 



2 The enthalpy (free energy) of the cavities is estimated from the total AGN power as H = -^jPV, 
where the PV work done by the jet as it inflates the cavity is determined by measuring the cavity size and 
its surrounding pressure, and the adiabatic index 7 is related to the internal composition of the cavities, still 
unknown. In the case of a relativistic plasma 7 = 4/3 and H = APV. 
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diameters, suggesting they are in pressure balance with the surrounding medium. Since the 
work required to inflate the radio lobes generally exceeds the nonthermal minimum energy 
of the radio emission (e.g., De Young 2006 and references therein), in the hypothesis that 
the assumptions adopted into the equipartition calculation are correct, some additional com- 
ponent besides the relativistic plasma may fill the cavity and contribute to support them 
internally. We want to investigate what is currently providing the necessary pressure support 
to sustain the cavities. Finally, in a more general view, we want to investigate what is the 
potential impact that such powerful outbursts have on the global properties of the ICM. 
Understanding well the ICM physics is essential to use galaxy clusters as high-precision 
cosmological tools. 

We address these problems by studying the X-ray properties of the most energetic 
outburst known in a galaxy cluster. MS0735+7421 (hereafter MS0735) is at a redshift of 
0.216. With a Hubble constant of H = 70 km s _1 Mpc" 1 , and fi M = 1 - fi A = 0.3, the 
luminosity distance is 1069 Mpc and the angular scale is 3.5 kpc per arcsec. 

2. Observation and data preparation 

MS0735 was observed by XMM-Newton in April 2005 during revolution 0973. In this 
paper only data from the European Photon Imaging Camera (EPIC) are analyzed and dis- 
cussed, while those from the Reflection Grating Spectrometer (RGS) will be presented in 
a subsequent paper. The MOS and pn detectors were both operated in Full Frame Mode 
with the THIN filter, for an exposure time of 70.6 ks for MOS and 60.7 ks for pn. We use 
the SASv6.5.0 processing tasks emchain and epchain to generate calibrated event files from 
raw data. Throughout this analysis single pixel events for the pn data (PATTERN 0) are 
selected, while for the MOS data sets the PATTERNS 0-12 are used. The removal of bright 
pixels and hot columns is done in a conservative way applying the expression (FLAG==0). 
To reject the soft proton flares we accumulate the light curve in the [10-12] keV band for 
MOS and [12-14] keV band for pn, where the emission is dominated by the particle-induced 
background, and exclude all the intervals of exposure time having a count rate higher than 
a certain threshold value (the chosen threshold values are 0.4 cps for MOS and 0.7 cps for 
pn). The remaining exposure times after cleaning are 50.4 ks for MOS1, 49.5 ks for MOS2 
and 42.0 ks for pn. Starting from the output of the SAS detection source task, we make 
a visual selection on a wide energy band MOS & pn image of point sources in the field 
of view. Events from these regions are excluded directly from each event list. The source 
and background events are corrected for vignetting using the weighted method described in 
Arnaud et al. (2001), the weight coefficients being tabulated in the event list with the SAS 
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task evigweight. This allows us to use the on-axis response matrices and effective areas. 
Unless otherwise stated, the reported errors are at 90% confidence level. 

2.1. Background treatment 

The background estimates are obtained using a blank-sky observation consisting of 
several high- latitude pointings with sources removed (Lumb et al. 2002). The blank-sky 
background events are selected using the same selection criteria (such as PATTERN, FLAG, 
etc.), intensity filter (for flare rejection) and point source removal used for the observation 
events. This yields final exposure times for the blank fields of 365 ks for MOS1, 350 ks for 
MOS2 and 294 ks for pn. Since the cosmic ray induced background might change slightly 
with time, we compute the ratio of the total count rates in the high energy band ([10-12] 
keV for MOS and [12-14] keV for pn). The obtained normalization factors (1.266, 1.303, 
1.283 for MOS1, MOS2 and pn, respectively) are then used to renormalize the blank field 
data. The blank-sky background files are recast in order to have the same sky coordinates 
as MS0735. 

The usual approach to perform the background subtraction is described in full detail 
in Arnaud et al. (2002). This procedure consists of two steps. In a first step, for each 
product extracted from the observation event list an equivalent product is extracted from 
the corresponding blank- field file and then subtracted from it. This allows us to remove 
the particle background. However, if the background in the observation region is different 
from the average background in blank field data, this step could leave a residual background 
component. The residual background component is estimated by using blank- field-subtracted 
data in a region free of cluster emission (in particular we consider an annulus lying between 8 
and 10 arcmin) and then subtracted in a second step from each MOS and pn product. In our 
case the residual is negative in the energy band adopted for the spatial analysis ([0.4-2] keV, 
see Sect. 13.11) . The residual count rate summed over the three detectors is —0.07 counts/s, 
which represents ~ 10% of the total background count rate in this energy band. 

The fluorescent emission lineal excited by the energetic charged particles that pass 
through the detector exhibit spatial variation over the detector. This effect would com- 
promise the reliability of the second step of the background subtraction for spectra, as the 
resulting total background spectrum (sum of the particle background spectrum, estimated in 



3 Al and Si lines at 1.5 and 1.7 keV, respectively, in the MOS data, and Ni, Cu and Zn lines around 8 
keV in the pn data. 
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the blank- field-subtracted cluster region, and the residual background component, estimated 
in the blank-field-subtracted outer annular region) shows fluorescence lines which are shifted 
by a few energy channels with respect to the MS0735 spectrum. Therefore the residual 
background component is neglected for the purpose of the spectral analysis. 

3. Morphological analysis 

The adaptively smoothed, exposure corrected MOS1 count rate image in the [0.3-10] 
keV energy band is presented in Fig. [U The smoothed image is obtained from the raw 
image corrected for the exposure map (that accounts for spatial quantum efficiency, mirror 
vignetting and field of view) by running the task asmooth set to a desired signal-to-noise 
ratio of 20. Regions exposed with less than 10% of the total exposure are not considered. 

The inner part of the cluster shows the high surface brightness characteristic of a cool- 
ing flow. We notice a sharp central surface brightness peak at a position 07 h 41 m 44 s .06 + 
74°14 / 38".62 (J2000). At small radii, two strong depressions having a diameter of about 150 
kpc (see Sect. I3.ip are visible on opposite side of the cluster center in the NE-SW direction. 
The cavities are surrounded by a bright X-ray emission of elliptical shape, in the form ex- 
pected by a radio cocoon. The elliptical discontinuity in the X-ray surface brightness appears 
more evident in the Chandra data and has been interpreted as a weak shock (McNamara et 
al. 2005). Beyond the cavity region the cluster maintains a slightly elliptical morphology 
up to large radii, showing hints of structures that could be the results of past radio activity 
that created older cavities. 



3.1. Surface brightness profile 

We consider different sectors to the north and south (N sector: between 10° west and 
60° east of north; S sector: between 15° east and 60° west of south) and to the east and 
west (E sector: between 30° north and 75° south of east; W sector: between 30° south and 
80° north of west) in order to include and exclude the cavity regions, respectively. We also 
consider a full 360° sector in which we mask the cavities (hereafter undisturbed cluster). 
For each sector, we compute a background-subtracted, vignetting-corrected, radial surface 
brightness profile in the [0.4-2] keV energy band for each camera separately. For the pn 
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Fig. 1. — M0S1 image of MS0735 in the [0.3-10] keV energy band. The image is corrected 
for vignetting and exposure and is adaptively smoothed (signal-to- noise ratio = 20). 
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data, we generate a list of out-of-time eventsj (hereafter OoT) to be treated as an additional 
background component. The effect of OoT in the current observing mode (Full Frame) is 
6.3%. The OoT event list is processed in a similar way as done for the pn observation event 
file. The profiles for the three detectors are then added into a single profile, binned such 
that at least a signal-to-noise ratio of 3 is reached. 

The cluster emission is detected up to 1.3 Mpc (~ 6 arcmin). In Fig. [2] we show the 
X-ray surface brightness profile for the sector containing the northern cavity compared from 
that of the undisturbed cluster. We note that the data of the undisturbed cluster appear 
regular, while those for the N sector show a clear deficit of emission between radii ~ 50 and 
~ 200 kpc relative to the other directions. The S sector profile shows a behavior similar to 
the N one, although the depression is less pronounced, whereas the E and W sectors appear 
regular. 

The surface brightness profile of the undisturbed cluster is fitted in the CIAO tool Sherpa 
with various parametric models, which are convolved with the XMM point spread function 
(PSF). The overall PSF is obtained by adding the PSF of each camera (Ghizzardi 2001), 
estimated at an energy of 1.5 keV and weighted by the respective cluster count rate in the 
[0.4-2] keV energy band. A single /3-model (Cavaliere & Fusco Femiano 1976) is not a good 
description of the entire profile: a fit to the outer regions shows a strong excess in the center 
as compared to the model (see Fig. [3]). The centrally peaked emission is a strong indication 
of a cooling flow in this cluster. We find that for 100 kpc <r <1150 kpc the data can be 
described (xLi ~ 1-65 for 87 d.o.f.) by a /3-model with a core radius r c = 195 ± 4 kpc and a 
slope parameter j3 = 0.77 ± 0.01 (3 a confidence level). The single /3-model functional form 
is a convenient representation of the gas density profile in the outer regions, which is used 
as a tracer for the potential. This best fit model is thus used in the following to estimate 
the cluster gas and total mass profiles (see Sect. ED- 

We also consider a double isothermal /3-model and find that it can account for the entire 
profile, when the very inner and outer regions are excluded: for 10 kpc <r <1150 kpc the best 
fit parameters are r cl = 200 ± 4 kpc, (3 1 = 0.77 ± 0.01, r c2 = 156 ± 8 kpc, f3 2 = 4.67 ± 0.56; 
Xred ~ 1-87 f° r 96 d.o.f. By assuming a common (3 value we find: r cl = 215 ± 4 kpc, 
r c2 = 30 ± 2 kpc, = 0.79 ± 0.01; x? ed ~ 2.08 for 97 d.o.f. 



Out-of-time events are caused by photons which arrive while the CCD is being read out, and are visible 
in an uncorrected image as a bright streak smeared out in RAWY. 
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Fig. 2. — Background subtracted, azimuthally-averaged radial surface brightness profile of 
the N sector data in the [0.4-2] keV range. The dotted red line shows the profile in full 
360° sector with masked cavities (undisturbed cluster). A deficit of emission in the N sector 
between radii of ~ 50-200 kpc is visible. 
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Fig. 3. — Background subtracted, azimuthally-averaged radial surface brightness profile in 
the [0.4-2] keV range. The best fit /3-model fitted over the 100 - 1150 kpc region is over- 
plotted as a solid red line. In the central region, the observed surface brightness profile shows 
a strong excess as compared to the extrapolation of the model. 
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4. Temperature map 

The temperature image of the cluster central region shown in Fig. H] is built from X- 
ray colors. Specifically, we extract mosaiced MOS images in four different energy bands 
([0.3-1.0] keV, [1.0-2.0] keV, [2.0-4.5] keV and [4.5-8.0] keV), subtract the background and 
divide the resulting images by the exposure maps. A temperature in each pixel of the map 
is obtained by fitting values in each pixel of these images with a thermal plasma, fixing Nh 
to the Galactic value (Dickey & Lockman 1990) and the element abundance to 0.4 solar (see 
Sect. 15.11) . Besides the evidence that the very central region is cooler than the surrounding 
medium, we do not notice any particular structure in the temperature distribution. 

The regularity of the temperature distribution points to a relaxed dynamical state of 
the cluster, thus excluding the presence of an ongoing merger. Since cluster merging can 
cause strong deviations from the assumption of an equilibrium configuration, this allows us 
to derive a good estimate of the cluster mass (see Sect. ED- 



5. Spectral analysis 

Throughout the analysis, a single spectrum is extracted for each region of interest and 
is then regrouped to give at least 25 counts in each bin. The data are modeled using the 
XSPEC code, version 11.3.0. Unless otherwise stated, the relative normalizations of the 
MOS and pn spectra are left free when fitted simultaneously. We use the following response 
matrices: ml_534_im_pall_vl .2.rmf (MOS1), m2_534_im_pall_vl .2.rmf (MOS2), 
epn_ff20_sY9_v6.7.rmf (pn). 



5.1. Global spectrum 

For each instrument, a global spectrum is extracted from all events lying within 6 
arcmin from the cluster emission peak, which corresponds to the outermost radius determined 
from the morphological analysis (Sect. 13. ip . We test in detail the consistency between 
the three cameras by fitting separately these spectra with an absorbed mekal model with 
the redshift fixed at z=0.216 and the absorbing column fixed at the galactic value (Nr = 
3.49 x 10 20 atom cm 2 , Dickey & Lockman 1990) and studying the effect of imposing various 
high and low-energy cutoffs. We find good agreement between the three cameras in the [0.4- 
10.0] keV energy range (kT = 4.64+°;^ keV for MOS1, 4.41+°;^ keV for MOS2, 4.35^ 
keV for pn). 
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Fig. 4. — Temperature map obtained by using 4 X-ray colors ([0.3-1.0], [1.0-2.0], [2.0-4.5], 
[4.5-8.0] keV) and estimating the expected count rate with XSPEC for a thermal mekal 
model, with fixed absorbing column N-r = 3.49 x 10 20 atom cm -2 and metallicity Z = 
OAZq. The metallicity gradient in the central region (see Fig. [6]) does not affect the map 
significantly. The map has a pixel size of 1.1 arcsec and is smoothed with a Gaussian of 
a = 16.5 arcsec. Superposed are the X-ray contours. The features outside the last contours 
are not significant, as they are mainly due to noise fluctuations. 
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The combined MOS+pn global temperature (in keV) and metallicity (as a fraction of 
the solar value, Anders & Grevesse 1989) derived from the best fit (xLi ~ 1-18 for 1947 
d.o.f.) are respectively: kT = 4A3^q^I keV, Z = 0.351o;q3Z q . The unabsorbed luminosities 
in this model (estimated from the average of the fluxes measured by the three cameras 
after fixing Nr = 0) in the X-ray ([2.0-10.0] keV) and bolometric band are respectively: 
Lx = 4.61±0.04 x 10 44 erg s -1 , = 1.05±0.01 x 10 45 erg s -1 , where the errors are given 
as half the difference between the maximum and the minimum value. 



5.2. Projected radial profiles: temperature and metallicity 

We produce projected radial temperature and metallicity profiles by extracting spec- 
tra in circular annuli centered on the peak of the X-ray emission. The annular regions are 
detailed in Table [H The data from the three cameras are modelled simultaneously using 
a simple, single-temperature model (mekal plasma emission code in XSPEC) with the ab- 
sorbing column density fixed at the nominal Galactic value. The free parameters in this 
model are the temperature kT, metallicity Z (measured relative to the solar values, with the 
various elements assumed to be present in their solar ratios, Anders & Grevesse 1989) and 
normalization (emission measure). The best-fitting parameter values and 90% confidence 
levels derived from the fits to the annular spectra are summarized in Table [TJ 

The projected temperature profile determined with this model is shown in Fig. [5j It 
shows a rise from a mean value of ~ 3.7 keV within 70 kpc to ~ 5.4 keV over the 150-630 
kpc region, then it declines to a value ~ 2.8 keV in the outskirts of the cluster (up to 1.3 
Mpc). 

The metallicity profile is shown in Fig. [6j a gradient is visible towards the central region, 
the metallicity increasing from Z ~ 0.27 Z Q over the 150-630 kpc region to Z ~ 0.62 Z inside 
the central 70 kpc. Due to the poor photon statistics that does not allow us to derive accurate 
measurements, we exclude the last two bins. 

We also perform the spectral fitting by leaving the absorbing column density as a free 
parameter and find little variation in the results. 

5.3. Deprojection analysis: temperature, density, pressure and entropy 

profiles 

Because of projection effects, the spectral properties at any point in the cluster are the 
emission-weighted superposition of radiation originating at all points along the line of sight 
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Fig. 5. — Projected X-ray gas temperature profile measured in the [0.4-10.0] keV energy 
range. 
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Fig. 6. — Projected X-ray gas metallicity profile measured in the [0.4-10.0] keV energy range. 
The last two bins have been excluded due to poor photon statistics. 
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through the cluster. To correct for the effect of this projection, we perform a deprojection 
analysis on the same annular spectra used in Sect. 15.21 by adopting the XSPEC pro jet 
model. Under the assumption of ellipsoidal (in our specific case, spherical) shells of emission, 
this model calculates the geometric weighting factor according to which the emission is 
redistributed amongst the projected annuli. 

The deprojection analysis is performed separately for the MOS1+MOS2 and pn spectra. 
The results are reported in Table 2 and the corresponding deprojected temperature and 
metallicity profiles are shown in Figs. [7]and[SJ respectively. We also perform the deprojection 
analysis by fitting simultaneously the spectra of the three cameras. In general, the results 
obtained from pn spectra appear more reliable, as the deprojected MOS and combined 
MOS+pn temperature profiles show some signs of instability and we do not manage to 
derive a temperature estimate in the last annulus. In the following general discussion and 
in the estimate of the cluster mass derived from the density and temperature profiles (see 
Sect. [6]) we therefore adopt the pn deprojection results. As expected, the deprojected central 
temperature is lower than the projected one, since in the projected fits the spectrum of the 
central annulus is contaminated by hotter emission along the line of sight. 

In Figs. [9lfTTlwe show various quantities derived from the deprojected spectral fits. The 
electron density n e (Fig. [9]) is obtained from the estimate of the Emission Integral EI = 
J n e n p dV given by the mekal normalization: 10~ 14 -E , //(47t[Da(1 + z)] 2 ). We assume n p = 
0.82n e in the ionized intra-cluster plasma. By starting from the density and temperature 
information derived from the deprojection analysis, we can calculate the pressure profile 
(Fig. [TO!) as P = nkT, where we assume n = 2n e . The average pressure surrounding the 
cavities is ~ 6 x 10~ n erg cm -3 . Similarly, the entropy profile (Fig. [TTI) is calculated from the 
temperature and density profiles by using the commonly adopted definition S = kTric 2 ^ 3 . 
Allowing for the different resolutions of Chandra and XMM, the high central entropy level 
of about 40 keV cm 2 that we measure is in agreement with the value of about 30 keV cm 2 
measured with Chandra (McNamara et al. 2005). As pointed out by Voit et al. (2005), this 
is consistent with being the result of the accumulation of the very energetic kinetic power 
outburst at the cluster center. 

We investigate the effect of changing the radial binning by performing a similar spectral 
analysis (projected and deprojected) on different annular regions. We find results consistent 
with those presented in Sections 15.21 and 15.31 We also note that the results presented here 
are in agreement with those derived from the analysis of Chandra data (McNamara et al. 
2005). 
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Fig. 7. — Deprojected X-ray gas temperature profile measured from MOS1+MOS2 (blue 
triangles) and pn (red full circles) spectra. For comparison, the projected profile is also 
shown (black open circles). 
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Fig. 8.— Deprojected X-ray gas metallicity profile measured from MOS1+MOS2 (blue 
triangles) and pn (red full circles) spectra. For comparison, the projected profile is also 
shown (black open circles). The last two bins have been excluded due to poor photon 
statistics. 
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Fig. 9. — Electron density profile measured from pn spectra (full circles). For comparison, 
the profile measured from MOS1+MOS2 spectra is also shown (open triangles). The solid 
line indicates the density profile obtained from the /9-parameters derived by fitting the surface 
brightness profile over the 50 - 1000 kpc. See text for details. 
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Fig. 10. — X-ray gas pressure profile measured from pn spectra (full circles). For comparison, 
the profile measured from MOS1+MOS2 spectra is also shown (open triangles). The red 
triangle represents the pressure of the thermal plasma filling the northern cavity, as estimated 
in Sect. 15.51 
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Fig. 11. — X-ray gas entropy profile measured from pn spectra (full circles). For comparison, 
the profile measured from MOS1+MOS2 spectra is also shown (open triangles). 
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5.4. Cooling Flow 

MS0735 was identified as a candidate "cooling flow" cluster because of extended Ha 
emission from its central galaxy (Donahue, Stocke, & Gioia 1992) and ROSAT HRI de- 
tection of a central high surface brightness peak (Donahue & Stocke 1995). However, this 
identification still has to be confirmed by a detailed X-ray spectral analysis. 

The surface brightness profile, temperature map and temperature profile derived from 
our XMM data analysis give indications of the presence of a cooling core. Here we further 
investigate the physical properties of the ICM in the central region. The cooling time is 
calculated as the time taken for the gas to radiate its enthalpy per unit volume H using the 
instantaneous cooling rate at any temperature: 

t ~ H - 7 kT (1) 

n c nuA(T) 7 - 1 fiX n n e A(T) 

where: 7 = 5/3 is the adiabatic index; \i pa 0.61 (for a fully-ionized plasma) is the molecular 
weight; Ah ~ 0.71 is the hydrogen mass fraction; and A(T) is the cooling function. We 
calculate the electron density by following the procedure described in Sect. I6.2[ using the 
/3-parameters derived by fitting the surface brightness profile over the 50 - 1000 kpc region^ 
(data in this region can be approximated by a f3 model with r c ~ 150 kpc and a slope 
parameter (3 ~ 0.70). We note that in the radial range where the (3 model is a good 
representation of the observed surface brightness profile, the density profile derived from 
this method agrees with that derived from the deprojection analysis of spectra extracted in 
concentric annuli (see Fig. [9]). 

Following Birzan et al. (2004), we define the cooling radius as the radius within which 
the gas has a cooling time less than 7.7 x 10 9 yr, the look-back time to z — 1 for our adopted 
cosmology. With this definition, we find r coo i ~ 100 kpc which corresponds to the central 
30 arcsec. The accumulated spectrum within this radius is extracted and compared to three 
different spectral models. Model A is the mekal model already used in Sect. 15. 11 Model B 
includes a single temperature component plus an isobaric multi-phase component (mekal + 
mkcf low in XSPEC), where the minimum temperature, kT\ ow , and the normalization of the 
multi-phase component, Normi ow = M, are additional free parameters. This model differs 
from the standard cooling flow model as the minimum temperature is not set to zero. Finally, 
in model C the constant pressure cooling flow is replaced by a second isothermal emission 
component (mekal + mekal in XSPEC). As for model B, this model has 2 additional free 



5 The best fit obtained in Sect. 13.11 cannot be extrapolated to the central region (r <100 kpc) therefore 
cannot be used here for the purpose of calculating the central cooling time. 
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parameters with respect to model A: the temperature, kT\ ow , and the normalization, Normi ow , 
of the second component. 

The results, summarized in Table El show that the statistical improvements obtained 
by introducing an additional emission component (models B or C) compared to the single- 
temperature model (model A) are significant at more than the 99% level according to the 
F-test. With our data, however, we cannot distinguish between the two multi-phase models. 
This means that the extra emission component can be equally well modelled either as a 
cooling flow or a second isothermal emission component. We note that the fit with the 
modified cooling flow model sets tight constraints on the existence of a "pedestal" minimum 
temperature (~ 1.5 keV). The nominal mass deposition rate in this empirical model is 
~ 260 ± 30 M Q yr _1 . We also attempt the fit with a classical cooling flow model (model D), 
by imposing kT\ ow = 0.1 keV, and find a mass deposition rate M ~ 40 ± 10 M yr _1 . 



5.5. Cavity regions 

Several states of matter have been proposed to fill the cavities (Pfrommer et al. 2005), 
including a population of relativistic electrons radiating at low radio frequencies and a di- 
lute, shock-heated thermal gas. The synchrotron emission from radio-filled cavities provides 
evidence of the existence of relativistic electrons and magnetic fields. The cavities may also 
be filled with a shock-heated thermal gas that can contribute to their internal pressure. 
This possibility is strengthened if the assumption of equipartition holds. Indeed, since the 
equipartition estimates of the nonthermal pressure in the radio bubbles give values which 
are typically a factor of ten smaller than the thermal pressures of the surrounding X-ray gas 
(e.g., Blanton et al. 2001, De Young 2006), the fact that the cavities are long-lived indicates 
that the necessary pressure support might be supplied by an additional thermal component. 
Observationally, besides the detection of hot X-ray emitting gas claimed by Mazzotta et al. 
(2002) within the ghost cavity of the MKW3s cluster, there has been no detection of hot gas 
over the cavities even in the clusters with better data available. 

In the case of MS0735, the existence of a relativistic plasma filling the cavities is clearly 
indicated by the presence of the radio source at the position coincident with the holes in 
the X-ray emission (McNamara et al. 2005). In order to investigate the possibility of an 
additional thermal gas component in the cavities, we perform a detailed analysis by modeling 
the spectra extracted in the cavity regions as the sum of ambient cluster emission and a hot 
thermal plasma, each with a characteristic temperature. The northern and southern cavity 
regions, detailed in Fig. [12], contain 10083 and 8872 source counts with a count rate of 0.23 
and 0.20 counts/s (MOS+pn), respectively. The results of the spectral analysis are reported 
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in Table H 

In the northern cavity, which shows a higher brightness contrast than the southern 
one in the XMM image, we find indications of the presence of a ~ 13 keV plasma. By 
assuming that the emission of this hottest component comes entirely from the cavity region 
whereas that of the coolest component is due to the projected foreground and background 
cluster emission^, we can estimate the density of the thermal plasma in the cavity to be 
n e ~ 6.4 x 10~ 3 cm -3 . This leads to an estimate of its pressure of ~ 2.7 x 10~ 10 erg cm -3 , 
which is a factor ~ 3-4 higher than the pressure of the surrounding medium (see Fig. [TO]) . 
However, we note that the improvement in the fit due to the extra thermal component is 
only marginally significant^. We therefore cannot place strong constraints on the existence 
of a hot thermal component filling the cavities and the estimate presented above has to be 
considered with caution. 



5.6. Shock front 

Chandra observation of MS0735 reveals a discontinuity in the X-ray surface brightness 
that has been interpreted as a weak cocoon shock driven by the expansion of the radio lobes 
that inflate the cavities. The shock front has an elliptical shape, being located ~ 240 kpc 
and ~ 360 kpc from the cluster center in east-west and north-south directions, respectively 
(McNamara et al. 2005). We find indications of a surface brightness feature at a position 
consistent with that seen in the Chandra data, however we cannot clearly distinguish a 
break in surface brightness that would be a signature of the shock front. By considering the 
elliptical shock front position derived from Chandra data, we perform a spectral analysis on 
the regions detailed in Figs. [T3] and O in order to derive the temperature of the pre-shock 
and post-shock gas. The results are reported in Table [51 

The particular choice of the elliptical regions next to specific sectors of the shock front 
(see Fig. \TS§ allows us to compare directly the spectral properties of the pre-shock and 
post-shock regions in different directions with respect to the azimuthally-averaged profiles. 
We do not find any appreciable variation in temperature (see Fig. [151) . In the analysis 
of the Chandra data, the shock properties are determined using a spherical hydrodynamic 
model of a point explosion at the center of an initially isothermal, hydrostatic atmosphere. 
When adapted to the XMM data by taking into account the smearing effect due to the 



6 This assumption is justified by the fact that a component as hot as 13 keV is detected nowhere else in 
the cluster and therefore is plausibly located only in the cavity region. 

7 X 2 /dof = 385/311 compared to x 2 / dof = 394/313, with F-test probability=0.028. 
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Fig. 12. — Regions considered for the spectral analysis of the cavities. The N and S cavities 
have a radius of 105 kpc and are located at a distance of ~ 170 and ~ 180 kpc from the 
cluster center, respectively. 
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large PSF, such a model predicts a temperature rise of ~ 10% (see Fig. [16]). Our spectral 
results obtained from elliptical annular regions next to the whole shock front (see Fig. Hlj) 
do suggest a temperature jump at the shock front, in agreement with those derived from 
previous Chandra observations, although due to the large error bars the pre-shock and post- 
shock regions are still consistent with being isothermal (see Fig. [T71) . From this analysis 
we therefore conclude that the XMM data can neither confirm or deny the existence of the 
shock front. 



In the following we estimate the total mass of the cluster using the usual assumptions of 
hydrostatic equilibrium and spherical symmetry. Under these assumptions, the gravitational 
mass M to t of a galaxy cluster can be written as: 



where G and m p are the gravitational constant and proton mass and /j, ~ 0.61. The mass 
contribution from galaxies is small, thus we neglect it. Therefore M tot (< r) = M gas (< 
r) + Mdm(< r )i i-e. the total gravitational mass within a sphere of radius r is given by gas 
plus dark matter mass. 

The mass profile derived from the equilibrium equation is strongly dependent upon the 
measured temperature profile. Large errors and irregular radial distribution of the temper- 
ature values induce large scatter on the reconstructed gravitational mass measurements. In 
order to obtain reliable mass estimates it is therefore crucial to select bins with a robust 
temperature estimate. For this reason we exclude the last temperature bin, although its 
inclusion would allow us to trace the temperature profile and therefore mass profile at much 
larger radii. 

The deprojected d In n g /d In r is calculated from the parameters of the /9-model derived 
in Sect. 13.11 In particular, the advantage of using a /3-model to parameterize the observed 
surface brightness is that gas density and total mass profiles can be recovered analytically 
and expressed by simple formulae: 



6. Mass profile 



6.1. Total gravitational mass 




(2) 




(3) 
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Fig. 13. — Elliptical regions considered for the spectral analysis of the pre-shock and post- 
shock gas. The centers of the ellipses are chosen to be at the same radial distance from the 
cluster center (~ 300 kpc), in order to allow us to make a direct comparison of the spectral 
properties in different directions with respect to the azimuthally- averaged profiles (see Figs. 
ED and [27]). 
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Fig. 14. — Elliptical regions considered for the spectral analysis of the pre-shock and post- 
shock gas. The concentric annular regions are defined in order to surround the whole shock 
front as derived from the Chandra data. 
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Fig. 15. — The open circles represent the projected temperature profile measured in Sect. 
15.21 The blue triangles indicate the measurements in the pre-shock regions along the east- 
west direction, the red square indicate the measurements in the post-shock regions along the 
north-south direction (see Fig. [T3|) . The error bars are at la level. 



-30- 



MS0735. 6 + 7421 model k162_015 




UJ 



00 

d 



20 40 60 80 100 

r (pixel) 

Fig. 16. — Temperature jump expected by the shock model. The dotted line shows the 
projected temperature profile as it would be seen by Chandra, and the full line is the result 
after smearing with the XMM PSF. The temperatures are in units of the pre-shock temper- 
ature, assumed to be constant. The radial distance is in Chandra pixel unit (1 pixel = 0.492 
arcsec) . 
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Fig. 17. — Temperature measured in the annular pre-shock (blue) and post-shock (red) 
regions (see Fig. [141) . The dashed line represents the position of the shock front along the 
east-west direction. The error bars are at la level. 
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M (< r) - fcr2 [ 3f3rT - —] (4) 
Gfirrip r 2 + r 2 dr 

As a first-order approximation, the temperature gradient is estimated by dividing the de- 
projected temperature profile (see Fig. 0) into two radial intervals and least-squares fitting 
a straight line in each interval. The total gravitating mass distribution derived from Eq. H] 
is shown in Fig. [18] as a solid line, with errors coming from uncertainties in the temper- 
ature measurement and /5-model parameters represented as long dashed lines. In Fig. [18] 
we also show (dotted line) the total mass calculated by assuming a constant temperature 
of < Tx >= 4.79 keV, where the average emission- weighted cluster temperature < Tx > 
is derived by extracting the global spectrum of the cluster after excluding the cooling flow 
region (i.e. in the annular region 30 arcsec - 6 arcmin). We note that the total integrated 
mass within a particular volume is dependent upon the local physical properties (local tem- 
perature and density gradients) and is is not strongly affected by the regions interior, or 
exterior, to that radius. The mass profile derived with this method is thus reliable in the 
region where the /3-model is a good representation of the observed surface brightness profile 
(100 kpc < r <1150 kpc, see Sect. 13.11) . whereas it cannot be extrapolated to the central 
region. 

As an alternative method, we then calculate the total mass by making direct use of the 
deprojected gas temperature and electron density values estimated from the spectral best-fit 
with a single phase model (Sect. 15. 3[) . In particular, the total mass enclosed within the 
midpoint between two consecutive shells is calculated from the pressure gradient and the gas 
density at the midpoint by following the method described in Voigt & Fabian (2006). The 
mass profile reconstructed by using this method, is shown as full red circles in Fig. [181 The 
uncertainties are calculated using error propagation from the pressure and density estimates 
(90% confidence level). We note that the cavities may affect the mass determination. In 
particular, since the cluster does not exhibit a smooth pressure profile between 100-400 kpc 
(see Fig. [TO]) , which would result in a non-monotonically increasing mass profile in this 
region, we exclude the fourth data point in the pressure and gas density profiles. 

The mass measurements reconstructed through the two different methods (direct ap- 
plication of the hydrostatic equilibrium equation with density profile derived from either 
/3-model fit to the surface brightness profile or deprojection analysis on annular spectra) are 
in very good agreement (see also Fig. [211 . 
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Fig. 18. — Mass profiles calculated with different methods. Red circles: total mass profile 
derived from deprojected gas temperature and density values. Black solid line: total mass 
calculated from Eq. H]by estimating the temperature gradient from the deprojected temper- 
ature profile (with error on the mass calculation coming from the temperature measurement 
and /9-model shown by the black long dashed lines). Black dotted line: total mass calculated 
from Eq. H]by assuming a constant temperature of 4.79 keV. Red triangles: gas mass profile 
derived from deprojected analysis. Black dashed line : gas mass profile derived from /9-model 
fit. See text for details. 



6.2. Gas mass and gas mass fraction 



We estimate the gas mass through two different methods, i.e. by using the density 
profile derived from either the /9-model fit to the surface brightness profile or the deprojection 
analysis on annular spectra. To convert from electron number density to gas density we use 
the relation p = i.83/im H ?7. e . 

In the first case, the gas mass profile is derived by integrating the gas density given by 
Eq. [3] in spherical shells and using the /3-model parameters determined in Sect. 13.11 The 
normalization of Eq. [3] is obtained from the combination of the best-fit results from the 
spectral and imaging analyses, which allows us to determine the conversion count rate - flux 
used to derive the bremsstrahlung emissivity that is then integrated along the line-of-sight 
and compared with the central surface brightness value. The adoption of the parameters of 
the /9-model fit in the outer regions produces an underestimate of the derived central electron 
density (see Fig. ^ and this turns into an underestimate of the gas mass in the central shells. 
Since the integrated gas mass is calculated by summing from the center outwards, any error 
in the measurement at small radii will propagate out to larger radii. However, the gas mass 
at small radii is much less than at large radii and any uncertainty in the measurements in 
the core are unlikely to have a significant effect on the gas mass profile further out. The gas 
mass profile derived from the /3- mo del is shown as a dashed line in Fig. [TSJ The gas mass 
profile obtained by integrating in spherical shells the gas density derived directly from the 
deprojected electron number density is shown in Fig. [18] as red triangles. Again, we note 
that the two profiles are in very good agreement. 

The gas mass fraction is the ratio of the total gas mass to the total gravitating mass 
within a fixed volume: 



The gas mass fraction profile derived from M tot and M gas measured from the /9-model and 
deprojection methods are plotted in Fig. [19] as black lines and red points, respectively. 
Adding the mass contribution from galaxies (~l-2%, Lin et al. 2003) to the total mass 
would have a small effect on the estimate of / gas , with variations lying within the error bars. 



Simple assumptions for the formation and evolution of galaxy clusters predict self- 
similar scaling relations. Galaxy clusters are assumed to form by spherical collapse of dark 
matter and gas due to gravitational instability. The gas is heated by shocks and adiabatic 
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Fig. 19. — /g as profile derived with different methods. Red circles: / gas profile derived from 
M to t and M gas measured directly from the deprojected density and temperature values (see 
red circles and red triangles in Fig. US]) . Black solid line: / gas profile (with error shown 
by the black long dashed lines) derived from M tot and M gas measured by using the density 
profile derived from the /3-model fit to the surface brightness profile (see black solid line and 
black dashed line in Fig. CO 
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compression. For dissipationless collapse, the X-ray emitting gas can be considered to share 
the same properties as the dark matter (e.g. Bryan & Norman 1998; Arnaud & Evrard 1999). 
Self-similar dark matter haloes are assumed to have identical properties when scaled by the 
virial radius, which is the radius separating the region where the cluster is in hydrostatic 
equilibrium from where matter is still infalling. However, deviations from self-similarity are 
expected in the ICM due to the details of individual collapse histories and to any physical 
processes beyond simple gravity and gas dynamics (e.g. Evrard & Henry 1991; Bryan & 
Norman 1998; Borgani et al. 2002 and references therein). In particular, the ICM of cooling 
flow clusters is affected significantly by heating that boosts the thermal energy of the gas. 

Here we investigate the correlations between the physical quantities observed for MS0735 
in order to scale it in the overall populations of clusters and assess its conformity to the 
generally observed mass-temperature (M-T) and luminosity-temperature (L-T) relations. 
Our final aim is to evaluate the impact of energetic AGN outbursts on such scaling relations, 
which are the foundation to construct the cluster mass function and to use these virialized 
objects as cosmological probes. 



7.1. Modeling the observed mass profile: determination of r& 

In order to rank MS0735 in the overall population of clusters we investigate the relations 
among different physical quantities considering their values at a given overdensity A. This 
is defined with respect to the critical density p CtZ = (3H%) / (8irG) , and within a cluster 
described as a sphere with radius r&: 

a = 3Mt ° t(< : a) (6) 

where the Hubble constant at redshift z is equal to (e.g., Bryan & Norman 1998): 



H z = H 0V /a rn (l + zf + l-n m = H h{z) (7) 

In order to estimate for various overdensities we need to fit the mass profile with an 
analytical function, that allows us to extrapolate the mass profile beyond the outer radius 
accessible to our X-ray observations. As a cosmologically motivated dark matter mass model, 
we consider the integrated NFW (Navarro et al. 1996) dark matter profile: 

M DM « r) = i*rl^- ^ 1 + eW(1 + e / (8) 
The scale radius r s and the concentration parameter c are the free parameters. 
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We perform the fit to minimize the y 2 of the comparison between the mass predicted by 
Eq. [8] and the mass profile reconstructed from both the density profiles estimated from the j3- 
model and deprojection methods. The best fit parameters for the mass profile reconstructed 
from density values estimated from the /3-model are: r s = 498+gg kpc, c = 3.45^Q^g (% 2 /dof = 
0.08/4). The best fit parameters for the mass profile reconstructed from deprojected density 
values are: r s = 1343^240 kpc, c = 1.66io!2i (x 2 /dof = 0.18/3). The quoted error are at the 
68% confidence levels (ler). The comparison of the two best fits of the NFW mass profile 
is shown in Fig. The corresponding gas mass fraction estimated from the NFW mass 
profile fit is shown in Fig. [21] for both the /9-model and deprojection methods. Note that we 
neglect the gas mass contribution to the total mass and we assume M tot (< r) = M DM (< r). 
We also perform the same fitting procedure by including the gas mass, i.e. by assuming 
M tot (< r) = M DM (< r) + M gas (< r), and find very similar results. 

From the best fit parameters we compute ta from Eq. [6] for various overdensities: 
A = 2500,1000,500,200. For tt m + Q A = 1, the virial radius corresponds to r A if A = 
A vir = 178fi° n 45 (Lacey & Cole 1993; Eke et al. 1996, 1998). For Q m = 0.3 adopted here, 
A v i r = 104. The virial radius derived in an Einstein-de Sitter Universe (f2 m = 1) is still often 
used in the literature, corresponding to a density contrast of A vir = 178 (or A vir = 200), 
so that r vir ~ r 20 o- The relation r vir = c x r s holds for the NFW mass profile. The results 
obtained from the NFW fit to the mass estimated through the /3-model and deprojection 
methods are shown in the Tables [6] and [7j respectively. We also report the total mass and 
the gass mas fraction. The error related to the mass estimate is obtained from half the 
difference between the maximum and the minimum value calculated at each radius for the 
set of parameters acceptable at la. The gas mass fraction is calculated at given overdensities 
by adopting the gas mass computed from the integration of the /3-model fit to the density 
profile extrapolated out to r&. The reported error is obtained from the error propagation 
and is dominated by the error on the total mass. 

We note that the value obtained for r 2 oo is highly dependent on the accuracy of r s , since 
it involves extrapolating the NFW model using r s . The measure of r s is very uncertain as 
the best-fit scale radius lies beyond the radius at which mass measurement can be made 
and is therefore strongly dependent on the outermost datapoint of the mass profile. The 
value for r 2 oo is shown for interest, but is not used for any subsequent analysis. We instead 
adopt r25oo> which is the most reliable estimate since it lies well within the radius at which 
mass measurements are obtained. We indeed find an excellent agreement between the values 
of 

^2500, ^2500 and / gas ,2500 obtained from the NFW fit to the mass estimated through the 
/3-model and deprojection methods. In particular, in the following discussion we adopt the 
values obtained with the /3-model method as their determination is more precise (cf. Tables 
E and ED. 
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The gas mass fraction profile scaled in radial units of r25oo is shown in Fig. [221 The 
apparent flattening of / gas outside r 2 soo is in agreement with previously observed profiles 
(Allen et al. 2002, Voigt & Fabian 2006), although it is systematically slightly above them. 
In particular, we find /gas, 2500 = 0.165 ±0.040, which is higher than the average value derived 
in a number of previous measurements with Chandra (e.g., Allen et al. 2002, Vikhlinin et al. 
2006). The high central gas fraction is close to the global baryon fraction in the Universe, 
constrained by CMB observations to be Q b /Q m = 0.175 ± 0.023 (Readhead et al. 2004, 
Spergel et al. 2003). However, our estimate of the central / gas could be affected by several 
effects. In particular, the high measured central gas mass fraction could result from an 
underestimation of the total mass, supporting the presence of non-thermal pressure in the 
core. It could also result from an overestimate of the gas density due to the boost in emissivity 
produced by the cavities (see Sect. 17.61) . 



7.2. Scaled temperature profile 

In order to make a fair comparison between the physical properties and investigate the 
scaling relations it is important to correct for the effects of the central cooling flow when 
measuring the characteristic temperature of the cluster. The average emission-weighted 
cluster temperature is calculated by fitting with a mekal model the spectrum extracted up 
to the outer radius detected by our X-ray observation (6 arcmin), after excising the cooling 
region (30 arcsec). We find a value < Tx >= 4.79 ± 0.13 keV. 

We show in Fig. [23] the projected temperature profile of MS0735 (red triangles) mea- 
sured in the region internal to r 250 o- The profile is scaled at an overdensity of 2500 and 
overlaid to the scaled temperature profiles of a sample of 6 relaxed clusters observed with 
Chandra (Allen et al. 2001). In Fig. [2l]we present the same kind of comparison with a sam- 
ple of 12 relaxed clusters for which the temperature profile has been measured farther out, 
up to more than 2 r 250 o, revealing a clear general temperature decline in the outer regions 
(Vikhlinin et al. 2005). In this case the temperatures are scaled to the emission- weighted 
cluster temperature, excluding the central 70 kpc region usually affected by radiative cool- 
ing. To be fully consistent with the scaled profiles presented by Vikhlinin et al. (2005), 
we thus estimate the emission-weighted temperature used for the scaling by extracting the 
global spectrum of MS0735 in the annular region 20 arcsec - 6 arcmin. We find a value 
{kT = 4.71 ± 0.11 keV) slightly lower than the "true" average emission-weighted temper- 
ature, as expected since the effect of cooling flow is not completely corrected (the cooling 
radius in MS0735 is ~ 100 kpc). The overall temperature profile measured for MS0735 is con- 
sistent within the scatter of the profiles observed in relaxed cluster, although for r > 0.5 r 2 5oo 
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Fig. 20. — Best fits of the NFW mass profile calculated with different methods. Red dashed- 
dotted line: best-fit NFW to the total mass profile measured from the deprojected density 
and temperature values, with la error shown by the red dotted lines. Black solid line: best-fit 
NFW to the total mass profile measured from the density profile derived from the /3-model 
fit to the surface brightness profile, with la error shown by the black long dashed lines. 
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Fig. 21. — / gas profile estimated from the NFW total mass profile derived with different 
methods (the gas mass is estimated from the /5-model fit to the surface brightness profile). 
The line colors and styles are the same as in Fig. [201 
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Fig. 22. — / gas profile with the radial axis scaled in units of r 2 5oo- The line colors and styles 
are the same as in Fig. [201 
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its values tend to lie at the top of the distribution. 

These results indicate that the energetic outburst in MS0735 does not cause dramatic 
instantaneous departures from the average properties of the ICM, as it has not had a large 
impact on the large-scale temperature profile. Of the roughly 30 clusters from the Chandra 
archive showing evidence of AGN outbursts (Rafferty et al. 2006), three have outburts of 
comparable energy, MS0735 (McNamara et al. 2005), Hercules A (Nulsen et al. 2005a) and 
Hydra A (Nulsen et al. 2005b; Wise et al. 2006). These three large outbursts all have ages of 
~ 10 8 yrs. Dunn et al. (2005) find that ~ 70% of cooling flow clusters currently show signs 
of outbursts, implying that outbursts are active most of the time. If the large outbursts were 
confined to only these few clusters, their incidence and ages imply that the clusters have 
undergone tens of outbursts since they were formed. In that case, the total energy added 
to the ICM outside the cooling region should have had a marked effect, contradicting our 
findings for MS0735. Thus, the relatively high incidence of large outbursts is more likely to 
be due to them occurring ~ 10% of the time in a significant proportion of all cooling flow 
clusters, rather than occurring most of the time in ~ 10% of clusters. 



7.3. Scaled metallicity profile 

We studied the metallicity profile in order to investigate the possible effect of uplifting 
and outflows on the metal distribution. Fig. [25] shows the projected metallicity profile of 
9 cooling flow clusters observed with BeppoSAX, rescaled to the virial radius (De Grandi 
& Molendi 2001). A strong enhancement in the abundance is found in the central regions. 
Overlaid is the metallicity profile that we measure for MS0735 (red triangles). We note that 
it is fully consistent with the strong central enhancement in the abundance exhibited by 
cooling flow clusters, indicating that the outburst has not smoothed the metal gradient. 

How to interpret this result in the context of current models is unclear. Simulations of 
metal evolution due to mixing and uplifting by buoyant bubbles indicate that low-power, 
subsonic flows have a relatively weak impact on the metallicity gradients found in clusters 
(Briiggen 2002, Omma et al. 2004), whereas powerful outbursts with short active time can 
flatten the central metallicity gradients if the metals are not replenished quickly (Heath et al. 
2006). The average jet power of the outburst in MS0735, as estimated from the shock model, 
is 1.7 x 10 46 erg s -1 (McNamara et al. 2005) and the observed steep central metallicity profile 
is consistent with the findings of Heath et al. (2006) for jet simulations of comparable power. 
However, the cavities in MS0735 would increase the jet power significantly, in which case 
the dredge up by the rising bubbles is expected to have a significant effect in flattening the 
metallicity gradient. If such outbursts are rare and we assume no continuous central metal 
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2500 



Fig. 23.— Temperature profile measured for MS0735 (red triangles) overlaid onto the 
temperature profile observed for a sample of 6 relaxed clusters presented by Allen et al. 
(2001). The profiles for all clusters are projected and scaled in units of T2500 and r 25 oo- ^2500 
is estimated within r 2 5oo an d results T 250 o = 4.65 ± 0.08 keV for MS0735. 
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Fig. 24.— Temperature profile measured for MS0735 (red triangles) overlaid onto the 
temperature profile observed for a sample of 12 relaxed clusters presented by Vikhlinin et 
al. (2005). The profiles for all clusters are projected and scaled in radial units of r 2 soo- 

The 

temperatures are scaled to the cluster emission-weighted temperature excluding the central 
70 kpc regions. 
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injection, the fact that we observe a normal gradient could thus be interpreted as a lack of 
strong evidence for significant metal uplifting or mixing by the cavities in MS0735. If instead 
as discussed above all clusters go through a similar phase, the observed agreement would not 
be surprising as the large outbursts would contribute to create the global cluster properties. 
Note however that the existing models do not allow us to exclude the possibility of some 
dredging occurring in MS0735. Indeed, the predictions of numerical simulations depend 
strongly on the details of the jet (like opening angle, initial Mach number, duration), which 
can significantly affect the extent of mixing. A direct comparison with the metallicity profile 
observed in MS0735 would therefore require a tailoring of such simulations to the particular 
conditions of the outburst in this cluster. Furthermore, observations of broader metallicity 
profiles compared to the cD light profiles in cD clusters (De Grandi & Molendi 2001, David 
et al. 2001, Rebusco et al. 2005) indicate that some outward diffusion of enriched gas due 
to AGN or mergers is occurring. Powerful outburst as the one in MS0735 are more likely 
to produce such an effect. What we conclude from our analysis is therefore that the single 
outburst in MS0735 has not appreciably affected the metal profile compared to other objects, 
but we cannot exclude the possibility of some mixing. 

On the other hand, we note that when dealing with azimuthally-averaged profiles we 
cannot study the local distribution of metals. In order to further investigate the metal- 
licity distribution we extract the spectra in the regions indicated in Fig. [26] and measure 
the abundance by fitting them with a mekal model. The results are reported in Table [8] 
and shown in Fig. [27J where we also plot for comparison the azimuthally-averaged profile. 
Simulations involving uplift by the cavities show that material is entrained at the upper (or 
lower) surface of the rising bubble and in its wake (e.g., Churazov et al. 2001), and the 
resulting metal distribution is very elongated along the direction of the bubbles (Roediger 
et al. 2006). We note that the possibility of making an accurate comparison between the 
pattern of abundance that we observe and the prediction of these models is complicated by 
the poor photon statistics, which do not allow us to place definite constraints on the metal 
distribution. Detailed abundance maps of the central cluster regions are required for this 
purpose. 



7.4. r - T relation 

A tight size-temperature relationship Is expected, as both of these quantities reflect the 
depth of the gravitational potential well in a virialized halo. Self-similarity predicts that 
ta oc y/<Tx> (Mohr & Evrard 1997). Studies of ensembles of simulated clusters have 
confirmed this proportionality, providing a value for the normalization in the relation (e.g., 
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Fig. 25. — Metallicity profile measured for MS0735 (red triangles) overlaid onto the metal- 
licity profile observed for a sample of 9 cooling flow clusters presented by De Grandi & 
Molendi (2001). The profiles for all clusters are projected and scaled in radial units of r v ; r . 
The virial radius is estimated from the relation r v ; r = 3.95 Mpc yj < Tx > / 10 keV (Evrard 
et al. 1996; see also Sect. 17.41) . H = 50 km s _1 Mpc" 1 , fi = 1, A = is assumed. 
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Fig. 26. — Regions considered to study the abundance distribution (cf. also the regions in 
Figs. H2and[l3D. 
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Fig. 27. — The abundances measured in different regions are indicated with the same color 
scheme used to define the regions in Fig. [26j For comparison, the azimuthally-averaged 
metallicity profile is also shown (open circles). The error bars are at la level. 



Evrard et al. 1996). Recent XMM-Newton observations of a sample of six relaxed galaxy 
clusters hotter than 3.5 keV find for an overdensity A = 2500 that (Arnaud et al. 2005): 



( 



5keV 



) 



0.50±0.03 



h(z) r 2500 = (500 ±5)kpc 



(9) 



where the factor h(z) corrects for the evolution expected in the standard self-similar model 
(see Eq. [7j). In the case of MS0735 the scaling factor is small, being h(z) — 1.11 at z — 0.216. 
This relation allows us to estimate a value of r 2 5oo ~ 440 ± 5 kpc, which is consistent with 
the value of r 2 soo = 463 ± 65 kpc determined from the overdensity profile (see Sect. 17.11) . 
We therefore find no evidence that the outburst is visibly affecting the size-temperature 
relationship. 

We note that different studies on samples of galaxy clusters find that the r^-T relation, 
as deduced from simulated clusters (Evrard et al. 1996), deviates systematically from the 
measured values of r&, as inferred from the overdensity profile (Sanderson et al. 2003; 
Piffaretti et al. 2005). The largest discrepancy is found for A = 200 and in the smallest 
haloes. Our result, derived for A = 2500 in a massive cluster, does not contradict these 
findings. The estimate of r 2 5oo is more reliable than that of r 2 oo as it does not require much 
extrapolation of the universal NFW profile (see Sect. 17.11) . Furthermore, the impact of 
additional additional, non-gravitational heating is minimal in massive clusters. Its effect is 
most pronounced in poor, cool clusters, where the extra energy required to account for their 
observed properties is comparable to their thermal energy (e.g., Ponman et al. 1996; Tozzi 
& Norman 2001). 



Theoretical predictions based on the simplistic assumption of an isothermal sphere for 
both the gas represented by its temperature and the collisionless dark matter particles (e.g. 
Kaiser 1986; Bryan & Norman 1998) infer a self-similar scaling relation between M tot and 
T at a given overdensity in the form M oc T 3 / 2 . Various observational studies have found 
different and sometime conflicting results regarding the slope and normalization of the M-T 
relation (e.g., Allen et al. 2001; Finoguenov et al. 2001; Ettori et al. 2002; Sanderson et 
al. 2003, Arnaud et al. 2005 and references therein). The relation derived by Arnaud et 
al. (2005) for a sub-sample of six relaxed clusters hotter than 3.5 keV observed with XMM- 
Newton is consistent with the standard self-similar expectation, following the relation: 
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M-T relation 
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5keV 
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1.51±0.11 



h(z)M 2500 = (1.79 ± 0.06) x 10 14 M ( 
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This result is in agreement with Chandra observations (Allen et al. 2001). In the case of 
MS0735, Eq. [Hturns into an estimate of M 2500 = (1.51±0.05) x 1O 14 M . By considering the 
whole XMM-Newton sample (ten clusters in the temperature range [2-9] keV), the relation 
steepens with a slope ~ 1.70 (Arnaud et al. 2005) indicating a breaking of self-similarity. 
In this case we estimate M2500 = (1.41 ± 0.05) x 1O 14 M . Although in better agreement 
with the M-T relation predicted from the cluster scaling laws, the mass estimate of M2500 = 
(1.74±0.42) x 1O 14 M as derived from the overdensity profile (see Sect. 17. II) is still consistent 
with a steeper slope of the relation considering the large errors that we measure. 



7.6. L - T relation and pre-heating 

The observed relation between X-ray luminosity and gas temperature in clusters is 
steeper than expected if cluster growth were governed by gravity alone. This steepening is 
best explained by the addition of heat to the ICM, and it has been estimated that the excess 
energies required to reconcile models and observations lie in the range of 1 to 3 keV per 
particle (Wu et al. 2000). The additional non gravitational heating is thought to have been 
injected into the gas during an early epoch of star formation and AGN activity, and therefore 
it has been called "pre-heating". Despite this term, there are now strong indications that 
powerful AGN outbursts occurring at late times may contribute a significant fraction of the 
extra energy (McNamara et al. 2005; Nulsen et al. 2005a, 2005b). In particular, in the case 
of MS0735 the driving energy of the shock as determined using a spherical hydrodynamic 
model is E s pa 5.7 x 10 61 erg (McNamara et al. 2005). The AGN outburst in this cluster is 
therefore heating the gas mass within 1 Mpc (~ 7.7 x 1O 13 M ) at the level of about 1/4 keV 
per particle. The heating level increases to ~ 0.6 keV per particle when considering the gas 
mass within r 2 5oo (^f g as,25oo ~ 2.9 x 1O 13 M ). This is a substantial fraction of the 1-3 keV 
per particle of excess energy required to preheat the cluster. 

In order to place MS0735 in the observed luminosity vs. temperature relation for galaxy 
clusters we estimate that its X-ray luminosity in the 0.1-2.4 keV energy range is Lx ~ 
5.0 x 10 44 erg s . For a total, emission- weighted temperature of kT = 4.4 keV the observed 
luminosity is a factor ~ 2.6 higher than that predicted by Markevitch's relation (Markevitch 
1998). The departure of MS0735 from the mean L-T relation is reduced when we correct 
both temperature and luminosity for the effect of the cooling flow, in the usual manner for 
the study of the scaling relations. In order to be fully consistent with the results presented 
by Markevitch (1998), we estimate the effect of the cooling flow by adopting the same 
approach of this author. The corrected temperature is obtained by averaging the projected 
temperature profile after excluding the coolest component in the central 20 arcsec bin (as 
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estimated from a multiphase spectral fit), resulting kT ~ 5.1 keV. The corrected luminosity 
is estimated from the flux observed by masking the central 20 arcsec and then multiplying 
by a factor 1.06 which accounts for the flux of the hot plasma component inside the masked 
region. The corrected luminosity in the 0.1-2.4 keV energy range is Lx ~ 4.2 x 10 44 erg s" 1 
and the corrected bolometric luminosity is Lboi ~ 9.0 x 10 44 erg s _1 , which correspond to a 
factor ~ 2 higher than expected from the corrected average temperature of the cluster. 

Note that the cooling flow region in MS0735 is bigger (see Sect. 15. 4p than the average 
value adopted by Markevitch (1998) for the cluster sample, therefore the luminosities and 
temperature estimated above are not completely corrected for the effect of cooling flow. A 
more precise correction is obtained by estimating the luminosities from the spectra extracted 
after excluding the central 30 arcsec and then adding back in the luminosity expected from a 
/9-model profile inside the masked region. In this case we find values of Lp.1-2.4] ~ 3.8 x 10 44 
erg s _1 and Lboi = 8.1 x 10 44 erg s _1 . Similarly as done above, the temperature is estimated 
by averaging the temperature profile after excluding the coolest component in the central 
30 arcsec bin, resulting kT ~ 5.4 keV. By considering these values "fully corrected" for the 
effect of cooling flow, the departure from the observed L-T relation is slightly reduced to a 
factor ~ 1.7. The corrected luminosity vs. temperature relations for luminosities estimated 
in the X-ray and bolometric bands are plotted in Figs. [28] and [29], respectively. The MS0735 
representative points are shown as red full triangles when corrected consistently with the 
method adopted by Markevitch (1998), and as blue open triangles when fully corrected for 
the effect of cooling flow. The results on temperatures and luminosities estimated in this 
section are summarized in Table [9] 



7.6.1. The " -cavity effect" 

The energetic outburst and the consequently rising cavities uplift the central cool, low- 
entropy gas up to large radii, and the same time the compression in the shells increases the 
ICM density. This effect of cooling flow, cavities and bright shells (that for simplicity in 
the following we refer to as the 'cavity effect') results in an increase of emissivity and thus 
luminosity. By considering a simple phenomenological model of the gas emissivity which 
assumes that all the gas filling the cavities is compressed into the bright shells due to the 
cavity expansion, we estimate that the luminosity is boosted by a factor which depends 
upon the cavity radius and shell thickness. For the particular configuration of the cavities 
observed in MS0735, we expect an increase in luminosity by a factor as high as about 25% 
(see Fig. [30]), consistent with our measurements. Note that since the cavities lie outside 
the cooling region, their effect is not taken into account by the methods adopted above to 
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correct for the cooling flow. The emission from the cavity region should also be masked 
before estimating the average properties of the "undisturbed" cluster. We thus extract the 
spectra in the annular region 1.3 - 6 arcmin (i.e. we excise the central 1.3 arcmin region 
which includes the cavities). From these spectra we estimate: Lp.1-2.4] ~ 3.1 x 10 44 erg s _1 , 
Lhoi — 6.4 x 10 44 erg s _1 , where the missing luminosity expected from a /3-model profile 
inside the central masked region is added back in. A comparison of these values with those 
obtained by fully correcting for the effect of cooling flow allows us to evaluate the impact of 
the cavity effect on the luminosity (cf. blue open triangles and stars in Figs. |2"51 and l2T?j) . 
The observed luminosity increase is consistent with that expected by our simplified model 
of the structure and emissivity of the gas. 

We want to stress that the cavity effect may contribute to partially explain the upward 
departure of MS0735 from the mean L-T relation, accounting for <25% of the increased 
luminosity. It may also be relevant for the measurements of the gas mass fraction in galaxy 
clusters: the cavities create lumps in the ICM and this results in an increase of emissivity 
and therefore in an overestimate of the gas density. We find indeed some evidence that the 
gas mass fraction in MS0735 might be higher than typically observed (see Sect. I7.ip . 

8. Conclusions 

The main conclusions of this work can be summarized as follows: 

• The energetic outburst in MS0735 can heat the ICM up to a level of about 1/4 - 1/2 
keV per particle, depending on its radius of influence, thus contributing a substantial 
fraction of the 1 to 3 keV per particle of excess energy required to preheat the cluster 
(Sect. I7.6p . Only a few outbursts of this magnitude erupting over the life of a cluster 
would be required to preheat it. 

• MS0735 is a factor ~2 more luminous than expected from its average temperature on 
the basis of the observed L-T relation for galaxy clusters (Sect. 17. 6p . Flux limited 
samples of distant X-ray clusters may be biased in favor of detecting clusters with 
energetic outbursts. 

• The ICM compression in the bright shells due to the cavity expansion produces a boost 
of the central luminosity of a factor that we estimate should be <25%, consistent with 
our measurements (Sect. I7.6p . This 'cavity effect' may contribute to partially explain 
the upward departure of MS0735 from the mean L-T relation, and may also lead to 
an overestimate of the gas mass fraction in qualitative agreement with the high value 
of /gas,25oo that we measure for MS0735 (Sect. 17.11) . 
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Fig. 28. — Luminosities in the 0.1-2.4 keV band corrected for the effect of cooling flow 
in the central ~ 70 kpc vs. emission-weighted temperatures derived excluding cooling flow 
components, from Markevitch (1998). The red full triangle represents MS0735 data from 
the present observations corrected consistently with the method adopted by Markevitch 
(1998). The blue open triangle represents the MS0735 data fully corrected for the effect of 
cooling flow (the cooling flow region in MS0735 is bigger than the average value adopted by 
Markevitch (1998) for the cluster sample). We also show for interest the data corrected for 
the cavity effect (blue open star). See text for details. 
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Fig. 29. — Same as Fig. [28] for corrected bolometric luminosities. 
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Fig. 30. — Estimated luminosity boost factor due to cavity espansion and gas compression 
in the shells, as a function of the thickness of the shell. Different curves refer to different 
radius Ri of the cavities (units in kpc). In this simplified phenomenological model all the 
gas filling the cavities is assumed to be compressed into the bright shells. For an adiabatic 
index 7 = 5/3, the maximum compression in a normal shock is a factor / = 4, so that the 
ratio of the shell thickness to the cavity radius is R2/R1 — 1 > (1 + 1/ f) 1 ^ 3 — 1 — 0.077. For 
the cavities observed in MS0735 (R± ~ 100 kpc), the shell thickness is ~ 8 kpc, leading to a 
luminosity boost factor <25%. 
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• The observed scaled temperature and metallicity profiles are in general agreement with 
those observed in relaxed clusters (Sects. 17. 21 and IT731 . Also, the quantities we measure 
for MS0735 are consistent with the r-T and M-T relations predicted by the cluster 
scaling laws (Sects. 17.41 and 17751) . This probably indicates that violent outbursts such 
as the one in MS0735 do not cause gross instantaneous departures from the average 
cluster profiles and cluster scaling relations (other than the L-T relation). However, if 
they are relatively common they may play a role in shaping these global properties. 

It is nevertheless hard to draw general conclusions from only one object, and therefore 
in the future we plan to observe and study more supercavity systems. 
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Table 1. Spectral fitting in concentric annular regions 



Radius 


kT 


Z 


xVdof 


(kpc) 


(keV) 


(Z Q ) 





0-70 (0-20") 
70-140 (20"-40") 
140-210 (40"-l'.0) 
210-315 (1'-1'.5) 
315-420 (l'.5-2') 
420-630 (2'-3') 
630-945 (3'-4'.5) 
945-1260 (4'.5-6') 



Note. — Results of the spectral fitting in concentric 
annular regions in the [0.4-10.0] keV energy range ob- 
tained by fixing the absorbing column density to the 
Galactic value (JV H = 3.49 x 10 2 °cm- 2 ). The temper- 
ature (in keV) and metallicity (in fraction of the solar 
value, Anders & Grevesse 1989) are left as free pa- 
rameters. Error bars are at the 90% confidence levels 
on a single parameter of interest. 



3.68lg;ig 0.62+°;°; 694/641 

4.53+°-^ 0.41+^07 695/652 

5.4S^M 0.321°$ 715/617 

5.651211 0.27±g;g| 718/690 

B.Utoil 0.241°;°° 484/505 

5.321°;^ 0.25±g;l| 582/599 

3.83lg;|| 0.44l°; 2 ^ 858/739 

2.80l};^ l.Hlngl 1129/810 



Table 2. Deprojection analysis 







MOS1 + MOS2 








pn 




Radius 


kT 


Z 


norm 


n e 


kT 


z 


norm 




(kpc) 


(keV) 




(xl0~ 4 ) 


(cm -3 ) 


(keV) 


(z & ) 


(xl0~ 4 ) 


(cm -3 ) 


D-7D 
70-140 
140-210 
210-315 
315-420 
420-630 
630-945 
945-1260 


3 65 + S ?S 

o.vo_ 19 

o c 7 +0.33 
°-° ' -0.28 

6.2415;" 
6.10lg;g 
5.00lg;P 
6.00l ;°J 
3.56l°;i 


o.so±8:ig 

n 77+O.23 
u - ' ' -0.20 

0.211™ 

0.371™ 
21 +0 - 20 

o.3ol°:l 
o.soll 


7 q7+0.36 
' -°'-0.36 

6.6418;^ 
8.86181 7 
12.8718:1 
8-6518:53 

7 i 9 +0.49 
' - iz -0.48 

5 29 +0 - 61 

u -^ y -0.59 


0210 +0 ,-°,°,^ 

u.uz,±u„g 0046 

o.oo75i8:88^8 
0.005318:881^ 

0035+ - 0007 

U.UUOO_ .0007 

0.002018:888^ 
0.000918:888^ 
o.ooo4i8:888? 


3 20 +0 /;f 

°-^ u -0.18 

4 1O+0.43 
^• la -0.39 

4.861°;™ 
5.80+8;?? 

4.9oi8;l 

6.3711;^ 
4.8311;!° 

9 71 +1.70 
z - 1 L -1.01 


59 + nTo 

u - uy -0.12 
u -°°-0.14 

n 47+0- 23 

u -^' -0.22 

0.221™ 
o.25i8;l? 
o.i2l8:g 
o.i5l8:?I 
i.59lS;» 


8 45+°^ 

y - dJ -0.58 
7 crc+0.55 
' •'- , '- , -0.48 

ll.63l8;i 

8-6618:1 

6 23+ - 37 
u - zo -o.48 

5.3418J 5 

o.98l8:^ 


0225 +0 /°, ^ 

U.U^Z J_q 0046 

o.oo9oi8;88^ 
o.oo49i8;881^ 

noss+o- 0007 

U.UUOO_ .0007 

o.oo2oi8;888^ 
0.000918:888^ 

04+ 0001 

0001 +°- 0001 

u.uuui_ .oooi 



Note. - Results of the deprojection analysis on annular MOS and pn spectra using the XSPEC 
proj ctxtbabsxmekal model. The column density is fixed to the Galactic value and the normalizations are 
in units of 10- 14 n c n p V/An[D A (l + z)} 2 . The fits give x 2 /dof = 3160/2610 and 2786/2635 for MOS and pn, 
respectively. 
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Table 3. Cooling flow analysis 



Par. Mod. A Mod. B Mod. C Mod. D 



kT 


°-° ' -0.09 


7.6ii; : | 


6 07+ 1 - 27 

°- u '-0.63 


4-25±8:i| 


Z 


U -°'-0.05 


0.52±g;8i 


^2 +0 - 05 

u - oz -0.06 


0.59lg;« 


Norm 


1.7 x IO" 3 


7.9 x 10™ 5 


9.1 x 10~ 4 


1.4 x 10~ 3 


kT\ ow 




1 AK+0-23 


9 97+O.4I 
' -0.39 


0.1 


Norm low 




M = 260 


6.7 x 10~ 4 


M = 40 


X 2 /dof 


891/787 


839/785 


839/785 


861/786 



Note. - The best-fit parameter values and 90% confi- 
dence limits of the spectral analysis in the central 30 arc- 
sec region. Temperatures (kT) are in keV, metallicities (Z) 
as a fraction of the solar value and normalizations in units 
of 10- u n c n p V/A7v[D A (l + z)f as done in XSPEC (for the 
mkcf low model the normalization is parameterized in terms 
of the mass deposition rate M, in M yr _1 ). See text for 
details. 
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Table 4. Cavity analysis 





Par. mekal mekal+mekal 


Cavity N 


kTi 5.24±g;|| 13.1™ 
Z 0.3ll£n 0.28±g;i| 
Norni! 4.9 x 1(T 4 2.3 x 10~ 4 
kT 2 3.45±S;gg 
Norm 2 2.7 x 10~ 4 
X 2 /DOF 394/313 385/311 


Cavity S 


fcTi 5.001°-! 0.43lg-^ 
Z 0.28t°i? 0.28i°;{ 4 
Normi 3.9 x 10" 4 

/cT 2 5.31_ ; 36 
Norm 2 3.8 x 10" 4 
X 2 /DOF 288/287 275/285 



Note. — The best-fit parameter values and 90% 
confidence limits of the spectral analysis in the cav- 
ity region. Temperatures (kT) are in keV, metal- 
licities (Z) as a fraction of the solar value and nor- 
malizations in units of 10 _ 14 n e n p V / 4n[D a(1 + z)] 2 
as done in XSPEC. 



Table 5. Shock analysis 



Region 


source counts 


kT 


Z 


X 2 /dof 




(MOS+pn) 


(keV) 


(Z e ) 





PRs-E 


3081 


6 1^+ L25 

u - id -0.88 


f) f ?4+°-33 


143/122 


PRs-W 


2361 


4 82 +0 - 71 


n 44+0-36 


87/100 


POs-N 


2958 


5 29+ - 72 


n 97+O.26 
u - z ' -0.22 


116/110 


POs-S 


3122 


r p-y+0.94 

o.oi _ 74 


25 +0 - 26 

U.ZO_ 24 


113/118 


(PRs-E) + (PRs-W) 


5442 


K Q4+0.58 
d -°^-0.50 


n 42+0-23 

u -^ z -0.20 


301/227 


(POs-N) + (POs-S) 


6080 


K 40+O.6O 


25+ - 17 

U - Zd -0.16 


233/233 


PRs (annulus) 


13187 


K 94+0-36 


21 +tuo 

u - zi -0.09 


440/461 


POs (annulus) 


22655 


5 63 +0 - 26 

O.UO_ 24 


29+ ' 07 

u - zy -0.07 


647/634 



Note. - Results from the spectral fitting in the regions indicated in 
Figs. [13] and [HI The fit is performed in the [0.4-10.0] keV energy range 
by fixing the absorbing column density to the Galactic value. When fitting 
simultaneously two different regions at the same radial distance from the 
center (e.g., PRs-E + PRs-W), the normalizations of each camera are 
linked in order to have the same value in the two regions. Error bars are 
at the 90% confidence levels on a single parameter of interest. 
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Table 6. Results from /5-model analysis 



A 




M tot 


/gas 




(kpc) 


(1O 14 M ) 




200 


1718 (218) 


7.12 (1.96) 


0.188 (0.052) 


500 


1107 (145) 


4.77 (1.26) 


0.181 (0.048) 


1000 


775 (105) 


3.27 (0.83) 


0.176 (0.046) 


2500 


463 (65) 


1.74 (0.42) 


0.165 (0.040) 



Note. - Characteristic radii r&, total mass 
Mtot and mass gas fraction / gas for various over- 
densities A derived from the NFW fit to the mass 
estimated through the /3-model. All the quanti- 
ties are estimated within r& (la errors in paren- 
theses). 
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Table 7. Results from deprojection analysis 



A 




M tot 


/gas 




(kpc) 


(1O 14 M ) 




200 


2230 (650) 


15.6 (8.78) 


0.109 (0.061) 


500 


1340 (410) 


8.46 (4.42) 


0.124 (0.065) 


1000 


875 (280) 


4.71 (2.33) 


0.141 (0.070) 


2500 


465 (158) 


1.76 (0.81) 


0.165 (0.076) 



Note. — Same as Table [6] for various overden- 
sities A derived from the NFW fit to the mass 
estimated through the deprojection method. 
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Table 8. Measured abundances 





Region 




r / 

z 

(Zq) 


X /dof 




PRs-E 




o 34 + S™ 


143/122 




PRs-W 




n 44+S-36 

u -^-0.30 


87/100 

O 1 / -L\J\J 




POs-N 




n 9 7+0.26 
U.Z 1 _ 22 


1 1fi/1 10 




POs-S 




25+S-25 

U.z,j_q 24 


113/118 

1 / ilu 




Cav-N 




o.3i±8.-g 


394/313 




Cav-S 




0.28j£? 


288/287 




A 






188/174 




B 




0.3018:2 


173/210 




C 




n Qfi+0.09 

u - OD -o.io 


346/380 




D 




0.37^;n 


371/335 


(PRs- 


-E) + (PRs-W) 


u -^ z -0.20 


301/227 


(POs- 


-N) + (POs 


-s) 


n 25+ - 17 

U - Zd -0.16 


233/233 


(Cav- 


-N) + (Cav- 


■s) 


n ^n+ - 08 

u - ou -0.08 


869/605 




A + B 




40 +0 - 13 


497/389 




C + D 




f) Q fi +0.07 
u - oo -0.07 


902/720 



Note. — Abundance measured in the regions 
indicated in Fig. [261 The mekal fit is performed 
in the [0.4-10.0] keV energy range by fixing the 
absorbing column density to the Galactic value. 
When fitting simultaneously two different re- 
gions at the same radial distance from the cen- 
ter (e.g., C + D), the normalizations of each 
camera are linked in order to have the same 
value in the two regions. Error bars are at the 
90% confidence levels on a single parameter of 
interest . 



Table 9. Results for L-T relation 



Band 



Data 




M 
rcl 



(keV) 



erg/s) 



0.1-2.4 
0.1-2.4 
0.1-2.4 
0.1-2.4 
Bolom. 
Bolom. 
Bolom. 



total 

CF corr. M 
CF fully corr. 
cavity effect corr 
CF corr. M 
CF fully corr. 
cavity effect corr 



4.4 ±0.1 
5.1 ±0.4 
5.4 ±0.4 



5.1 ±0.4 
5.4 ±0.4 



5.00 
4.19 
3.76 
3.06 
8.98 
8.12 
6.37 



4.20 
4.69 



1.89 
2.07 
2.26 



Note. — Temperatures and luminosities measured for MS0735. The 
effect of the cooling flow (CF) has been evaluated by employing two 
different approaches. CF corr} 11 indicates the correction adopted by 
Markevitch (1998): the temperature is obtained by averaging the tem- 
perature profile after excluding the coolest component in the central 
20 arcsec bin (as estimated from a multiphase spectral fit); the lumi- 
nosities are estimated from the flux observed by masking the central 
20 arcsec and then multiplying by a factor 1.06 which accounts for the 
flux of the hot plasma component inside the masked region. CF fully 
corr. indicates a more precise correction which takes into account the 
"true" extent of the cooling region in MS0735: the temperature is esti- 
mated by averaging the temperature profile after excluding the coolest 
component in the central 30 arcsec bin; the luminosities are estimated 
from the spectra extracted after excluding the central 30 arcsec and 
then adding back in the luminosity expected from a /3-model profile 
inside the masked region. Lf^ indicates the luminosity expected from 
the luminosity vs. temperature relation (Markevitch 1998). The lu- 
minosities corrected for the cavity effect are also reported. See text 



(Sect. 17101) for details. 



